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1 Introduction 



More than 40 years have passed since the Korteweg-de Vries (KdV) equation 
was solved by Gardner et al. [Ij using the inverse scattering method based on 
its Lax pair [2] . The number of known integrable systems following the KdV 
equation, particularly partial differential equations (PDEs) in 1 + 1 space- 
time dimensions, has increased enormously, and various techniques to study 
them have been developed. Recently the center of researchers' interest has 
shifted from continuous PDEs to differential-difference or partial difference 
equations wherein at least one of the independent variables takes discrete 
values. A major problem in this trend is how to find a suitable difference 
analogue of a given differential equation. The suitable discretization of an 
integrable continuous system is generally required to retain the integrabil- 
ity [3], but that is not sufficient if the original continuous system has some 
essential internal symmetries. This becomes conspicuous if we consider in- 
tegrable discretizations of nonlinear Schrodinger (NLS)-type systems, i.e., 
two-component systems of second order that allow the complex conjugation 
reduction between the two dependent variables. As a prototypical example, 
we discuss the cubic NLS system in 1 -|- 1 dimensions 



Note that the reduction r = aq* with a real constant a simplifies the two- 
component system (11. ip to the scalar NLS equation OE]. In addition, (11. ip 
is invariant under the space reflection x — )■ —x as well as the time reflection 
t ^ —t with q ■H- r. The suitable and elegant space discretization of the NLS 
system (II. ip was proposed by Ablowitz and Ladik [7J in the form 



Indeed, system (II. 2p is integrable and, with a rescaling of variables, reduces 
to (ll.ip in the continuous space limit. Moreover, (II. 2p allows the complex 
conjugation reduction between g„ and r„ and possesses the same invariance 
properties with respect to the space/time reflection as the continuous system 
(11.11) . Although 35 years have already passed since their work, the Ablowitz- 
Ladik discretization (II. 2p is still a rare example of success. Indeed, even now, 
only a small number of suitable space discretizations of integrable NLS-type 
systems are known (see, e.g., [8J); they are all (1 + l)-dimensional systems 
with only one discrete spatial variable. The problem of how to discretize the 



iqt + qxx - '^q^r = 0, 
in - r^^ + 2r^g = 0. 



(1.1a) 
(1.1b) 



iqn,t + (gn+1 + Qn-l " 2g„) - g„r„(gn+l + Qn-l) = 0, 

irn,t - (rn+i + r„_i - 2r„) + r„g„(?^n+i + ^n-i) = 0. 



(1.2a) 
(1.2b) 
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continuous time variable in such systems is an interesting topic [H], but we 
do not discuss it in this paper. 

The main objective of this paper is to provide the first example of a 
suitable discretization of an NLS-type system in 2 + 1 dimensions. In par- 
ticular, we consider the discretization of both spatial variables in a (2 + 1)- 
dimensional NLS system known as the Davey-Stewartson system [TD] (also 
see [H]). Note that the Davey-Stewartson system is integrable [T2HIS] and 
appears to be the only genuinely (2 + l)-dimensional generalization of the 
NLS system (11. ip (cf. the Calogero-Degasperis system [IZj). Moreover, even 
if we include other types of integrable systems, the list of known systems with 
two discrete and one continuous independent variables is still very short. 
Thus, it is a highly nontrivial and challenging task to obtain the suitable 
space discretization of the Davey-Stewartson system. To solve this problem, 
we first propose a natural (2 + l)-dimensional generalization of the Ablowitz- 
Ladik lattice (cf. (II. 2p ) by constructing its Lax pair. This (2 + l)-dimensional 
Ablowitz-Ladik lattice certainly reduces to the Davey-Stewartson system in 
the continuous space limit. A relevant Lax pair as well as the resulting sys- 
tem was previously studied by other authors [TH] (also see [H]), but the time 
part of our Lax pair is essentially more general than the previously known 
one [18]. As a result, the time evolution of our system is a linear com- 
bination of four elementary time evolutions, two of which were previously 
unknown. Moreover, it can be shown that the four time evolutions are mu- 
tually commutative. Thus, the (2 + l)-dimensional Ablowitz-Ladik lattice 
is general enough and appears to be promising. However, it does not allow 
the complex conjugation reduction directly and thus is not a suitable space 
discretization of the Davey-Stewartson system in its present form. To fix 
this shortcoming, we only have to consider a certain nonlocal transformation 
of dependent variables, which symmetrizes the equations of motion. Thus, 
we obtain the suitable space discretization of the Davey-Stewartson system 
that indeed allows the complex conjugation reduction between the new vari- 
ables after the transformation. In addition, the invariance properties of the 
continuous Davey-Stewartson system with respect to space/time reflections 
turn out to be properly incorporated in our space-discrete Davey-Stewartson 
system. 

This paper is organized as follows. In section 2, we propose a (2 + 1)- 
dimensional version of the Ablowitz-Ladik lattice by considering an appropri- 
ate generalization of the Lax pair for the original Ablowitz-Ladik lattice. To 
uncover how the complex conjugation reduction can be imposed as an NLS- 
type system, we consider a nonlocal change of dependent variables; it can 
turn the (2 -|- l)-dimensional Ablowitz-Ladik lattice into the suitable space 
discretization of the Davey-Stewartson system that indeed allows the com- 
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plex conjugation reduction between the new dependent variables. In section 
3, we elucidate how the general time evolution considered can be decomposed 
into four elementary time evolutions corresponding to the four directions on 
the two-dimensional lattice. On the basis of this decomposition and using the 
Hirota bilinear method [20], we construct some exact solutions of the (2 + 1)- 
dimensional Ablowitz-Ladik lattice and the space-discrete Davey-Stewartson 
system. In particular, multidromion solutions are presented explicitly. The 
last section, section 4, is devoted to concluding remarks. 



2 Derivation based on Lax pairs 

2.1 (2 + 1) -dimensional Ablowitz— Ladik lattice 

As a generalization of the Lax pair introduced by Ablowitz and Ladik [7], we 
consider the following linear system on the two-dimensional lattice: 

^jiV'i,™. • '^n+l,m '^n,m Qn,m(pn,m: (2.1a) 
^m^n,m ■ 0n,m+l 4^n,m '^n,m'^n,m) (2.1b) 
dipn m 

(2.1c) 
(2.1d) 

Here, and denote the forward difference operators in each spatial 
direction, and the parameters a and b are arbitrary constants. The time 
dependence of the functions is usually suppressed. The compatibility condi- 
tions (9tA+Vn,m = A+c^tV'n.m and (9t A+0„,m = A+(9t0„,m for the linear sys- 



5 



tern (12. ip provide the (2 + l)-dimensional Ablowitz-Ladik lattice, 



9(in,m r) _i_ h 

r\, ~r Qn+l,mJ-^n+l.m ~r "Q'n— 1 m (^Qn,m+1 ^n+l,mQn,m—l 

~l~ Qn,mBn,m -^n+l,mQn,m 0, (2.2a) 

Or^ m 

TT, ~l~ '^n,m+lCn,m+l (^'^n,m—l ^'"n+l,m Dn m+l'^n—l^m 

~\~ ^n,m^n,m -Srt,m+l'"n,m 0, (2.2b) 

^71+1, m ^n,m ^ijln,m+l^n,m Qn,m^n,m—l\ (2.2c) 

Bn.m+1 Bfi^m n+1 ,m(ln ,m ^n,mQn~l,m)^ (2. 2d) 

(1 Qn,rn^n,m)C'n^m Cfi^i^^ni^^ Qn,m—l^n,m—l)y (2.2e) 

(1 ^n,mQn,m)Dyi,fn Dn^n+liX ^n—l,mQn—l,m)' (2.2f) 

Indeed, if all the functions depend on n and m only through n + m, (12. 2p 
reduces to the original Ablowitz-Ladik lattice [?]; the latter contains the 
integrable discrete NLS system (11.21) as a special case. To restore the (1 + 1)- 
dimensional Lax pair involving the spectral parameter z, we set V'n.m = z'^'^ip'^^ 
and (f)n,m = z"^'"-'^ (j)'^ ^, rewrite the hnear problem (12.10 in terms of ^ and 
(pnmy then consider the dimensional reduction. 

Under appropriate boundary conditions at spatial infinity, we can use 
(I2.2cp - (]2.2fj) recursively to express the auxiliary fields An^m, Bn,m, Cn,m, and 
Dn^rn globally in terms of and r„ Thus, they can be considered as 
the defining relations for the auxiliary fields, and the (2 + l)-dimensional 
Ablowitz-Ladik lattice (12.21) has intrinsically nonlocal nonlinearity. Note that 
relations (I2.2cp - fl2.2fl) already appeared in the literature on an integrable time 
discretization of the (1 + l)-dimensional Ablowitz-Ladik lattice [3| [2T| [22]. 
Incidentally, in the stationary case of dtqn,m = dtfn,m = 0, (12. 2p reduces to a 
nontrivial system of partial difference equations in 1 + 1 dimensions. 

Using the simple transformation 

where 7 is a constant, and omitting the prime, we can introduce the terms 
+J(ln,m and —'jrn,m in (I2.2ap and (12.2b[) . respectively. Clearly, these terms 
can be absorbed by constant shifts of An^rn and Bn,m- 

2.2 Continuum limit 

By choosing the parameters appropriately and taking the continuous space 
limit, we can reduce system (12. 2p to the continuous Davey-Stewartson sys- 



6 



tern. To see this, we first shift the auxihary fields as 
and rewrite fl2.2l) as 

(2.3a) 

-|- fl^rn^rrt+l ~l~ ''"n,m— 1 '^^n,m) ^(^n+l/m ~l~ ^n— l,m 2r. 



dt 

(2.3b) 

-^n+l.m -^n.m ^{,Qn,m+l^n,in Q'n,m^n,m— l) ; (2.3c) 

Bn,m+1 Bn^m b{TnJ^\^rnQn,rn ^n,mQn—l,m)i (2.3d) 

(2.3e) 
(2.3f) 

Subsequently, we rescale the variables and parameters as 
qn,m = '^x ■ q{x,y), r^^m = Ay ■ r{x,y), x:=nAx, y -.= 171 Ay, 

^?i,m = ^(3^) y)) Bn^m = B(x,y), Cn^m = C {x , y) , Dn.m = D(x,y), 

wherein the time dependence is suppressed and 

b ^ 



{Ayf {Axf 

Thus, in the continuum limit Ax, Ay — >■ 0, ( 12. 3 p reduces to the continuous 
Davey-Stewartson system [I31II1IIS1I221I21] , 

qt + l^qxx - aqyy -{A + C)q + q{B + D) = 0, (2.4a) 

n + aryy - /3r,, + r{A + C) - (5 + D)r = 0, (2.4b) 

A^ = C^ = -a{qr)y, (2.4c) 

By = Dy = -/3(rg),. (2.4d) 



Note that the Davey-Stewartson system fl2.4p is a hnear combination of the 
two commuting flows corresponding to a = 0, /3 7^ and a 7^ 0, /3 = P^[2B] 
(also see [21] )• In subsection 13. H we present its discrete analogue, that is, 
the (2 + l)-dimensional Ablowitz-Ladik lattice fl2.2p is a linear combination 
of four commuting flows. This is a quite natural result because (i) each of 
the two Davey-Stewartson flows provides an asymmetric (2 + l)-dimensional 
generalization of the NLS system and (ii) the Ablowitz-Ladik discretization 
of the NLS system is actually a sum of two elementary flows (and one trivial 
flow) in the same hierarchy [31171 [27H29] . 



2.3 Noncommutative extension 

Actually, the (2 + l)-dimensional Ablowitz-Ladik lattice (12. 2p is integrable in 
the general case where the dependent variables take their values in matrices, 
as long as the operations such as addition and multiplication make sense. 
In that case, "1" in (I2.2ep and ( ]2.2fp should be interpreted as the identity 



matrix. 

We can further generalize it to a variable-coefficient system wherein the 
parameters a and h become arbitrary matrix- valued functions of one spatial 
variable as am '■= a{m) and 6„ := b{n). To obtain such an extension, we 
consider the following generalization of the linear system (12. ip : 



'4'n+l,m — ^'4'n,m + (ln,m4^n,my (2.5a) 

^ 4'n,m ~l~ ^n,m'^n,m,y (2.5b) 

~l~ ^ Qn~l,mbn~l'Pn,m Qn,mDn,rn'Pn~l,mj (2.5c) 

~l~ ^^n,m—l^m—l'^n,m ^n,mC'n,m'^n,m—l' (2.5d) 



d(f)n,m 

dt 



Note that the "spectral parameter" z is nonessential in the (2 + l)-dimensional 
case and can be fixed at 1 as described in subsection 12.11 The compatibility 
conditions for the linear system (12. 5p indeed provide the noncommutative 
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system with site-dependent coefficients, 

dc[n,m n _|_ h 

~l" Qn,mBn,m ^n+l,mQn,m ~ 0, (2.6a) 

~l~ '^n,m+lCn,m+l ~l~ ^n,7?i— l^m— 1 ^n'"n+l,m -^n.m+l'^n— l,m 
~l" ^n,m^n,m „j_|_ir„ „j = 0, (2.6b) 

^n+l,m ^71, m (^mQn,m+l^n,m ~\~ Qn,m^n,m—l(^m—ly (2.6c) 

Bn,m+1 Bfi fyi i^n^n+l,rnQTi,m ^n,mQn—l,rnbn—ly (2.6d) 

(-^ Qn,m^n,m)C'n^m C'^^-X ^(-^ Qn,in—l^n,in—l)i (2.6e) 
(-^ ^n,mQn,rn)DyiYn D^.m+li^I ^n—l,mQn—l,rn)- (2.6f) 

If and r„ ^ are rectangular matrices, the identity matrix / in (I2.6el) and 
that in (12.6f[) have unequal sizes. In the commutative case of the parameters, 
the site-dependent nature of (12.61) is nonessential if both nm=-oo ^^"^ 
n^-oo ^" nonzero finite values. Indeed, if we change the variables as 



m— 1 \ / 11— 1 \ / m— 1 \ / n— 1 



-1 



\]=-oo 



c — r n — n 

the site-dependent parameters can be normalized to 1. We can also obtain a 
similar result in the noncommutative case. 



2.4 Appropriate change of dependent variables 

For simplicity, in the following discussion, we consider only the commutative 
and constant-coefficient case wherein the parameters a and b are constants 
and all the quantities are scalar. Thus, the lowest-order conservation law for 
(12.21) is given by 

dlog{l-qn,mrn,m) _ r_, , D n _ ^"^ 1 

I ^Qn~l,m^n,m ~r J-^n,m\^ Qn—l,mf"n—l,m) Q.n,mTn—l,m\ 

~^ [ ^Qn,m^n,m—1 ~l~ C'n,m(l Qn,m—l^n,m—l} Qn,m—l1^n,m\ 

(2.7) 

The existence of an ultralocal conserved density log(l — qn,mrn,m) implies 
that a nonlocal transformation involving infinite products of (1 — qn,mrn,mY 
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with 5 7^ could be applied (cf. [30]); this is indeed the case as we will see 
below. 

The (2 + l)-dimensional Ablowitz-Ladik lattice (12. 2p is invariant under 
a space reflection (n, m) — )• (— m, — n) with a minor redefinition of the pa- 
rameters and the auxiliary fields. However, (12. 2p does not allow the complex 
conjugation reduction between g„ ^ and m in the local form. Therefore, 
we need to identify new "conjugate" variables instead of qn,m and r„ ,„ and 
rewrite (12.21) in a more symmetric form using the new variables. For this 
purpose, we consider a gauge transformation so that the spatial part of the 
Lax representation obtains invariance with respect to the combined space 
reflection (n, m) — )■ (— n, —m) or, equivalently, a 180 degree rotation around 
the origin. Thus, we apply the gauge transformation 



■i!) = X 6 =Y ^ 



{21 



to (I2.1al) and (12.1b|) and change the dependent variables as 



Ur. 



X~±IL^IIL1 ~ IL^IIL -17- 



(2.9) 



Here, X„ ^ and F„_m are defined as 



Xri 



^ n — 1 



j,mi 



Yr 



n,m 



m—1 



f n - (ln,krn,k. (2.10) 



The norming functions hm,{t) and /„(t) are introduced to realize the complex 

conjugation reduction between Un,m and Vn,m] they will be determined later. 

-1 ' 



One can also use {Y[t=ni \/l - qn,krn,k) ' instead of [IfcL-oo a/T" 
to maintain the invariance under the space refiection {n,m) — )■ (— m, — n). 
This modification causes no essential difference in the following discussion, 
so the transformed system can become ideally symmetric with respect to 
space reflections. Here and hereafter, we assume that \qn,mTn,m\ ^ 1 so that 
■\/l — <ln,mrn,m and its inverse as well as their inflnite products as consid- 
ered above are uniquely and well deflned. For example, we consider that 
qn,m = 0{Ax) and r^^^ = 0{Ay) (cf. (I2.1ap and fl2.1bp). and a/1 - qn,mrn,m 
is deflned as the Maclaurin series in qn,m.^n,m = 0{AxAy). Note that Un^mVn,m 



so that the inverse transformation of (12. 9p can be obtained immediately. 
Thus, the spatial part of the Lax representation acquires the form 



$n,m+l 



1 rn'^n,m 



yn^m. 1 



(2.11) 
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Very recently, D. Zakharov has considered essentially the same scattering 
problem in [SUES]- However, this is an accidental coincidence, because the 
first author arrived at this Lax representation as well as its generalization 
implied in subsection 12.51 independently before the papers [3T1[32] appeared. 
The invariance under the combined space reflection (n, m) — )■ (— n, — m) can 
be easily seen if we shift the indices of the linear wavefunction by 1/2, i.e.. 



n,m+ 



1 



1 



r 



where "ii^m =■ 'i'. 



and $r 



=: Indeed, because the determinant 



of the spatial Lax matrix above is unity, we obtain 



n- 



n,m+i 



It should be possible to apply the transformation (12. 9p with (I2.10p di- 
rectly to the (2 + l)-dimensional Ablowitz-Ladik lattice (12. 2p and derive the 
transformed equations of motion with the aid of the conservation law (12. 7p . 
However, the nonlocal nature of (12. 2p makes such a computation rather com- 
plicated and difficult. Thus, as an alternative, we apply the gauge transfor- 
mation (12. 8 p with (12. 9p and (I2.10p to the time part of the Lax representation, 
(I2.1cp and (]2.1dp , and determine the time evolution of the gauge-transformed 
wavefunction, dt'^n,m and dt^n,m- Its compatibility with the scattering prob- 
lem (12. lip can provide the transformed equations of motion. 



2.5 Space-discrete Davey— Stewartson system 

Before applying the transformation described in subsection 12. 4^ we fix the 
boundary conditions for the (2 -|- l)-dimensional Ablowitz-Ladik lattice (12. 2p 
as 



lim (g„,m,r„,m)= lim {qn,m,rn,m) = 0, (2.12a) 



m— )■— oo 



lim Cn,m = c ( ^?:!!^-J- ) , lim Dn,m = d ( ) • (2.12b) 

llm I m^-OO \ In / 



n— f — OO 



However, we do not fix lim„_s._oo An^m and limm^-oo -Bn,m in order to obtain 
interesting solutions such as dromion solutions; they can also depend on the 
remaining spatial variable and time t. In (I2.12ap . the dynamical variables 
qn^rn and rn,m are assumed to approach zero sufficiently rapidly. In (]2.12bp . 
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c and d are constants. The defining relations (I2.2ep and (12.2fl) enable C„ 
and Dn^m to be expressed globally as 



n,m— 1 / \ ^ n—l,'. 

Thus, the gauge transformation (12.81) with (12. 9p and (I2.10p changes (I2.1cp 
and fl2ldll to 



and 



d^n.m ■^n,m+l y _|_ 4 iT/ _l_ "i™ ,t, 

7TT ^ ^n,m+l ~r -'^n,m^ n,in ~r C-— ^n,m— 1 

Y Y 

, n,m I 7 n,m ^ 

~r "tT ^n—l,m^n,m ~\? ^n,m^n— l,m 

-'^n— l,m+l -'n— l,m 



rs, " -i^ ^n+l,m ^n,m^ n,im ~ \/- ^n—l,m 

^ n,m ^ n—l,m 

^n+l,m— 1 ^n,m— 1 



where 



-^n,m J-n,m 

Recalling that qn,m'rn,m = Un,mVn,m, the above relations combined with (12. lip 
comprise the Lax representation for the transformed system. To express it 
in a concise form, we introduce the quantities 

1 1 , , 

U!n,m ■= , = = , _ (2.14a) 



Jn,m '■— — — y-rn-l / . ' (/.i4DJ 

Yn+l,m _ nfc=-oo ~ ,kV n+1 ,k (0^A\ 

9n,m ■— yr — rrm-1 Aj ' (^/.i4Cj 
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Thus, we obtain the Lax representation in the form, 

^n+l,m ^n,m^n,m ~^ '^n,m'^n,m^n,my (2.15a) 
^n,m+l 'Wn,m^n,m ~l~ ^n,m^n,m^7i,m; (2.15b) 



dt 

dt 



(^fn,m,^n,m+l ~l~ ^n,m^ n,m ~l~ (^fn,m—l^n,m—l 

-\- bWn—l^m9n—l,m^n~l,m^n,m d(]n—l,m^n,m^n—l,my (2.15c) 

^fl'n,r?i^n+l,m ~l~ Bn^m^n,m ~l~ dcjn—l,m^n~l,m 
-\- (lWji ,iYi—lfn,m—l^n,rn—l^n,m ^fn,m—l^n,m^n,m—l- (2.15d) 

The corresponding boundary conditions are given by 

hm {Un,m,Vn,m) = lim {Un,m, Vn,m) = 0, 
n— ^— oo m— oo 

hm 1 . In 



hm fn^rn = 1 , lim gn,', 



7 ;yiL,m J 

ActuaUy, we can generahze fl2.15p to a more general form wherein the 
spatial part is given by 

^n+l,nJ ^n,m^ n,rn ~l~ Qn,m^n,mi 
^n,m-i-l Sji^m^n,m ~l~ ^n,'m^ n,mi 

with four independent functions Wn, qn,m, and rn,m- Thus, it is pos- 

sible to start with this general Lax representation and then consider the 
reduction. However, we skip such a discussion to maintain an easy-to-read 
flow of the paper. 

The compatibility conditions for the linear system (12.151) with Wn,m = (1 — Un,mVn,my 
provide the time evolution equations for Un,m and Vn^m- They can be written 
in a natural compact form using new auxiliary fields an,m and (3n,m defined 
as 

~ 1 1 

-^n,m • '^^n—l,m9n~l,m ip^n—l,m'^n,m diUn^rn^n—l,m) '^(-^n,rm (2.16a) 
Bn^m • '^Wn,m—lfn,m—l ('^^n,m^n,m— 1 CUn^rn—l1^n,m) '^(^n,m- (2.16b) 
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Thus, we finally arrive at the desired system, 

7 (f '^n,m'^n,rn) id'^n,m9n,in'^n+l,m ~l~ ^^n— l,mfi'n— l,m^7i— l,m 



dt 

CLWn^mfn,m'^n,m+l (^^n,m—lfn,m—l^n,m—l) 

(^n,m ~l~ '^n—l,m9n—l,m(fi^n—l,m^n,m ~l~ '^^n,m^n— l,m)] 0; 

(2.17a) 



b'Wn^m9n,m^n+l,m d'Wji—\^m9n—l,rnVn—l,rn) 



dt 

b'Wn^m9n,m^n+l,m d'Wji—\^m9n—l,rnVn—l,rn) 

1 

\ 

2 

Pn,m '^n—l,m9n~l.miP'^n—l,inVn,in ~l~ '^^ra,m'^n— l,m)] 0, 

(2.17b) 

,m,f 11,171 Wn,m+lfn+l,m if (o, c) 7^ 0, (2.17c) 

if (6,^)7^0, (2.17d) 

^n^'n-,m \^n,m—lfn,m—lifl'^n,m'^n,rn—l ~l~ (^'Un,rn—l'^n,'m)\ 1 (2.17e) 

)] • (2.17f) 

_ i_ 

Here, a, 6, c, and d are constants and Wn,m = (1 — Un,mVn,m) ^ ■ In the same 
way as (12.21) . (12.171) also admits a dimensional reduction to the Ablowitz- 
Ladik lattice [7j. Using f l2.17cp - fl2A7f|) . we can rewrite f l2.17ap and f l2.17bp 



in a more symmetric form with respect to space reflections. When c = —a* 
and d = —6*, the (2 + l)-dimensional system (I2.17P allows the complex con- 
jugation reduction Vn,m = o'u*^ ^ with a real constant a; in this reduction, 
the auxiliary fields fn,m and gn,m become real-valued, while the auxiliary 
fields an,m and /3n,m become purely imaginary. In particular, (I2.17P with 
purely imaginary a, b, c(= —a*), and d{= —b*) provides the suitable space 
discretization of the Davey-Stewartson system (cf. (12.41) ). 

Similarly to the continuous case (cf. [33H35] ). when c = —a and d = —b, 
we can consider the reduction of Vn,m = crun,m and a„,m = Pn,m = to obtain 
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a (2 + l)-dimensional analogue of the modified Volterra lattice [36j . 

dUn m 2 \ /? 7 

" \^ ^^nm) \^^n,m9n,mUn+l,m "^n— l,m5'n— l,m^n— l,m 



dt 

-\- OiWn^Ynfn,m'^n,m+l '^^n,m— 1 /n,m— l^n,m— l) ) (2.18a) 
,mfn,m 'Wn,m+lfn+l,m 

ifa^O, (2.18b) 

Wn,mgn,m = Wn+l,mgn,m+l if & 7^ 0. (2.18c) 

_ 1. 

Here, Wn,m = (l ~ '^'^n,m) ^ ■ would be interesting to look for a relationship 
between fl2.18p and the discrete modified Nizhnik-Veselov-Novikov hierarchy 
in 121] (also see gMH]). 

3 Solutions by the Hirota method 

In this section, we discuss how to construct exact solutions of the discrete 
Davey-Stewartson system (12.171) using the Hirota bilinear method [20]. Be- 
cause of the complexity and irrationality of the equations of motion, it would 
be too hard to solve fl2.17p directly, so we take an alternative approach. First, 
we bilinearize the (2 + l)-dimensional Ablowitz-Ladik lattice (12. 2p . Sub- 
sequently, we consider the effect of the nonlocal transformation (12. 9 p with 
fl230|) . f l213|l . fl2:ii|) . and fl2:T6|l in the bilinear formalism. The infinite prod- 
ucts appearing in the nonlocal transformation can essentially be expressed 
locally in terms of a "tau function" . Thus, we can obtain exact solutions of 
(12. 2p and (I2.17P concurrently from the same set of bilinear equations. 

3.1 Decomposition into four commutative flows 

Before applying the Hirota bilinear method, we demonstrate that the (2 + 1)- 
dimensional Ablowitz-Ladik lattice (12. 2 p can be decomposed into the four 
elementary flows. In view of (I2.2cp . fl2.2dp . and ( ]2.12bp . we rescale the aux- 
iliary fields as 



A =■ aA^°^ B =■ bB^°^ C =■ cC^ D =■ dD^^^ 



(3.1) 



The corresponding boundary conditions are 

lim {qn,m,rn,m)= lim {qn,m,rn,ni) = 0, 



m—^—OD 



lim = f^^V lim d(0) 

11111 ^ '7 I ' ^^^^^ ^n,m \ 7 

rim J m^-OD \ 
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Thus, considering the simplest cases where only one of the parameters a, 6, 
c, and d does not vanish, we obtain the four elementary systems: 



a-system 



dtaQn,m — (ln,m+l + ^n+l,mQn,m, (3.2a) 
dta^n.m ^n,m—l '"»i,m^rem' (3.2b) 

■^n+l,m ^n,m ijln,m+l'^n,m Q.n,mTn,m—l) i (3.2c) 



6-system 



dtfjQn,m — l,m ^n/m-S^J^j, (3.3a) 

dt^T^ .^ = T^^i^jn "I" -^i,rn+l'"n,m; (3.3b) 
■^n,m+l -^n)m ~ iTn+l,mQn,m '^n,m(ln-l,m) j (3.3c) 



c-system 



dt^qn,m = C*i+l,mQ'n,»Tt-l, (3.4a) 
dta^n,m = ~f'n,m+lCi,^m+ly (3.4b) 
(1 '~ (ln,m1^n,m)C^l^ = C^^i^j^{l — qn,m-l^n,m-l) , (3.4c) 



(i-system 



dt^qn,m — ~(ln+i,mDl^li^^, (3.5a) 

dt^1^n,m — -Oji^m+l'""— (3.5b) 
(1 ~ fn,m(ln,m)Dl^l^ = i5^j„_^]^(l — rn-l,m(ln-l,m) ■ (3.5c) 

Clearly, the time evolution in 02.21) is a linear combination of these four time 
evolutions, that is, dt = adt^ + bdt^ + cdt^ + ddt^. In fact, they are mutually 
commutative, so the above four systems belong to the same integrable hier- 
archy as the original system (12.21) . To check the commutativity conditions 
dudtpqn,m = dtpdt^qn,m and dt^dt^rn,m = dt^dt^rn,m for {«,/?} C {a,b,c,d}, 
we need to know how to express time derivatives of the auxiliary fields. Us- 
ing fl3.2l) - fl3.5l) . we can obtain all necessary expressions in the local forms, 
e.g., 

^ /l(0) _ r-)(0) _ 7-)(0) f 

^td n,m Hn,m+l-^n,m+l' n—l,m 'dn,'m-^n,m' n-l,in-li 
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Here, we assumed that all "integration constants" etc. can be set equal to 
zero. With these local expressions, we can check the commutativity of the 
four flows by direct computations. Note that lim„_^_oo dt^AnJn, hmm-s.-oo dt^B, 
lim„_> _ oo c?t^Cn°m, and \imm^_ao dt^DnJn do not vanish in general. Thus, the 
ta-fiow can change the boundary value of the auxiliary field in the tc-flow and 
vice versa; the same applies for the tfe-fiow and t^-fiow. 



3.2 Bilinearization 

Because the four systems (I3.2I) - (I3.5I) are compatible, in the sense that their 
flows mutually commute, we will consider here their common solution de- 
noted as qn,mita,tb,tc,td), rn,mita, t^, td) , ctc. Here, ta, h, tc, and td are 
independent arguments. Thus, the solution of the original system (12.21) is 
obtained by setting 

ta = at, tb = bt, tc = ct, td = dt, (3.6) 

which indeed implies the relation dt = adt^ + hdt^^ + cdt^ + ddt^ . 
We assume a solution expressible in the form, 

^n,m, Hn,m _ \ 

(ln,m ~ ~^ 5 fn,m — ~^ ; \O.I&) 

<L = log f = A+ (^-h^) , (3.7b) 

B^l = dt, log ( ^\ = K ( , (3.7c) 



n,m / \ n,m 



WW WW 
f~i{0) n,m+l-'^ n,m—l /")('-') ^ n+1 ,m-'^ n—1 ,m /q y i\ 



and bilinearize the four systems f l3.2p -f l33|) in terms of the "tau functions" 
Fn,m, Gn,m, and Hn,m as follows: 

• a-system 

Fn+l,m+l9taGn,m G'n,mf^ta ^n+l,m+l — Fn+l,mCrn,m+li (3.8a) 

Fn,mdta-^n,m Hn^^Ot^Fn^m Fn^rn+lHn,m—ly (3.8b) 

Fn,mdta-Pn+l,m Fn~\-\^mdtaPn,m GrnmHn,m—ly (3.8c) 

• 6-system 

Pn^mdtfjGji^^ Gn,mdti,F'n,m Pn+l,mGn—l,mi (3.9a) 

Fn+l,m+ldt^Hn,m ~ -f^n,m.'9fj,-^n+l,m+l = Fn,m+lHn+l,m, (3.9b) 

Fn,mdti,Ffi^m+l Fn^rn-\-ldtijFnYn G n—l^mFIn,m,^ (3.9c) 
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• c-system 

Fn+l,mdtfi n,m G n,mdtcFn+l,m — Fn+l,m+lG n,m—li (3.10a) 

Pn,m+l9tcFIn,m Hn^rn(^tcFn,m+l -^n,m-^n,m+l 5 (3.10b) 

Pn+l,mPn,m+l Pn+l,m+lPn,m CrnmHn^m: (3.10c) 

• d-system 

Pn+l,mdt^G fifii G fiYTidt^F^j^^fYi Pri,mGn+l,m^ (3.11&) 

Fn,m+ldt^H„ m ~ Hn ni^td-^n^m+l = -^n+l.m+l-^n-l.m; (3.11b) 

-^n+l,m-fn,m+l -fn+l,m+l-^n,m Gn^mHn^m- (3.11c) 

To be precise, each triplet of bilinear equations gives a sufficient condition for 
the corresponding original system. Note that for the c-system and (i-system, 
the bilinear forms as well as some exact solutions were studied in 1181. 



3.3 General solution formulas 

Once a solution of the bilinear equations fl3.8p - fl3.lip is obtained, formula 
fl3.7p with f l3.ip and fl3.6p provides the solution of the (2 + l)-dimensional 
Ablowitz-Ladik lattice f l2.2p . We assume that it satisfies the boundary con- 
ditions fl2.12p . Thus, by applying the nonlocal transformation (12. 9p with 
fl230|) . f l233|) . fl2Al|l . and fl^lH]) . we can also obtain the solution of the dis- 
crete Davey-Stewartson system fl2.17p . To evaluate the effect of this nonlocal 
transformation, we use fl3.7ap and fl3.10cp (or fl3.11cp ) to rewrite the infinite 
products as 



n— 1 n— 1 



nr. -TT / -^j+l,m+l-^j,m / -^n,m+l ^j,m 

V ^ j+l,m-'^ j,m+l y ^ n,m J ^ j,m+l 



m—1 m— 1 



nr. TT / -Tn+l./c+l-^n./c / -Tn+l.m ,. 



Fn+l,kFn,k+l y Fn^m k^-<^ -^n+l,A: 

Because we assumed \qn,mrn,m\ < 1, the value of {Fn+i,m+iFn,m) / {Fn+i-niFn,m+i] 
is always restricted to the neighborhood of 1. For simplicity, in considering 
the solution of (12.171) . we also assume that Fn^m is positive (or, at least, 
|argF„„i| is sufficiently small); the positivity condition > can fully 
justify the use of the formulas for the square root, such as y/X^ = X and 
^yX/Y = \/X/ \/Y. We set the norming functions and In in fl2.10p as 



hm = \ I lim _ , L = \ lim ^'^ 
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Thus, we obtain 



Y- _ I £n^m+l ^ _ I Fn+l,m ^ „n 

After all, we can express the transformation from (12.21) to fl2.17p locally in 
terms of the "tau function" Combining ([23]), f l27[3|) . (EUD, ([31]), f lHTj) . 

and (13.121) . we arrive at general solution formulas for the discrete Davey- 
Stewartson system (12.171) in the form. 



(3.13a) 



An,m = ^ {adt^ - bdt^ - cdt^ - rfS^J log ) , (3.13b) 



= - (69t, - adt^ - cdt, - ddtj log ' ) , (3.13c) 

\ ^ n..m 



/n,m = ; 5'n,m = ^ ■ (3.1od) 



Here, the time variables are set as in (13. 6p and the auxiliary fields an,m and 
f3n,m are determined from An^m and -Bn,m through (12.161) . Using the bilinear 
equations (I3.8cl) and (I3.9cl) and noting that (I2.17ep and (I2.17fj) are identities 



in a, b, c, and d, we obtain compact expressions for an,m and /3„,^ 



ttn.m = {~adt^ + cdtj log ( l^^p±l ) , (3.14a) 



Pn,m = i-bdt, + f/9,J log ( ) , (3.14b) 



-'■ n,r?i 

and new bilinear equations, 

Pn,mdtc-Fn+l,m Pn+l,mdtc-^n,in G^n,m~l-^n,m! (3.15) 
Fn,mdt^FYi^rn+l -^n,m+lf^tj-^n,m GnmHn—l,m- (3.16) 

Note that (I3.15P and (I3.16P fill in the piece missing in (I3.8p - (l3.1ip . In the next 
subsection, we construct common solutions to all these bilinear equations. 

As described below (12.171) . when c = —a* and d = —b*, we can impose 
the reduction Vn,m = o'u^ ^ with a real constant a. This reduction can be 
realized by requiring that Fn^m > and i?n,m = <^G*^^. 
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3.4 Solitons and dromions 



The set of bilinear equations fl3.8l) - fl3.11l) together with (13.151) and fl3.16p is 
not ideally symmetric in its present form. In particular, it is not clear why 
reductions such as F*^ = Fn,m and Hn,m = o'G* ^ are allowed. To restore 
the symmetry, we need only to rewrite fl3.8ap , (13.8bp , fl3.9al) , and (]3.9bp using 
dSScD, flTOj) . and flXnEj) (or fIXTO ). For example, using (IXScj) and then 
(I3.10cl) . (I3.8ap can be rewritten as 

Thus, the full set of bilinear equations can be reformulated in the symmetric 
form. 



Fn+l,mFn,m+l -^n+l,m+l-^n,m 


Gn.mHn^nii 


(3.17) 


■Fn,m+l(^taG n,m 


~ G n,m9ta-^n,m+l 


Pn,mGn,m+l} 


(3.18a) 


Fn+l,mdtaHn,m 




Fn+l,m+lFIn,m—li 


(3.18b) 


Fn,mdta-^n+l,m ~ 


Fn+l,mdta-^n,m 


G n,mH n,m—li 


(3.18c) 




"~ G yi ^ndtfjFYi m+l 


Fn+l,m+lG n—l,mi 


(3.19a) 


Fn+l,mdtfjFIji^m 


~ Hn,m9tijFj2-\-l^rn 


Fn,mHfi-i-irni 


(3.19b) 


F a F ^1 - 

n,m'^th-'- n,m+l 


- F ^1 a F = 


Gn~l,m,Hn,m} 


(3.19c) 


Fn+l,mdti, Gn^m 


Gn,m9tcFj2-\-l^rn ~ 


~ -^n+l,-m+lCn,m~l ; 


(3.20a) 








(3.20b) 


Pn,mdtc-^n+l,m 




Gn,m—lHn,mi 


(3.20c) 


Fn+l,mdt^G n^m 


G n,mdt^Fn-\-\^m ~ 


" Fji^mGn+l,mi 


(3.21a) 








(3.21b) 


Fn^mdt^F^jyij^Y 


Fn,m+ldt^F^,ryi 


Gn,mHn—l^m- 


(3.21c) 



It is now clear that the ta-flow and t^-flow can be identified with the tc- 
flow and t^-flow, respectively, through the complex conjugation reduction. 
Moreover, (13.180 and fl3.19p correspond to each other by the interchange of n 
and m, up to a redefinition of the variables. With these symmetries in mind, 
we can considerably reduce the task of constructing explicit solutions to the 
above 13 bilinear equations. 

In the same way as in the continuous case (see, e.g., |1Q]), we can construct 
the one-soliton solution and a two-soliton solution straightforwardly. The 
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one-soliton solution is given by 



{l-pp){l - qq) 

Cjn,m — gP Q e , -Hn,m — QP Q G , 

where cu := gt^ - + q'H^ - ptd, oJ := -q'Ha + ph - + P^^^d, and 
5f, p, g, etc. are nonzero constants. The constant q should not be confused 
with qn,m.- Actually, we can shift ta, ^6, ^c, and td in (13. 6p by arbitrary 
constants, but this freedom can be absorbed by rescaling g and 'g. When 
c = —a* and d = —b* (cf. (13.61) ). we set p = p*, q = q*, and 'g = erg* with 
(t(1 — |pp)(l — |gp) < 0. Thus, Fn^m > and -ff„,m = (^Cn m' complex 
conjugation reduction is realized. 

To save space, we omit a rather lengthy expression for a two-soliton solu- 
tion. These soliton solutions are direct (2 + l)-dimensional analogues of the 
soliton solutions of the (1 + l)-dimensional systems such as (II. 2p . With an 
appropriate choice of the parameters, they represent straight line solitons in 
the physical variables (cf. (I3.7al) or (I3.13al) ) and thus are not localized. 

In the following, we obtain more interesting solutions, that is, dromion 
solutions; dromions [H] are spatially localized "solitons" that decay expo- 
nentially in all directions |35] and can exhibit nontrivial interaction proper- 
ties [inill2]- More details as well as an extensive list of references can be 
found in the review article [13] . In analogy with the continuous case [ID] , the 
one-dromion solution is obtained as 

Fn,n. = 1 + W"e^^+^^ + ^ai^(gg)™e-^+^^ 

(1 -pp)(l - qq) 

^n,m — OLI2P Q e , rln,m — «2lP Q G , 

where 

ui := -p^^h-ptd, oJi := ph+p'^td, UJ2 := qta+q^h^, UJ^ ■= -q'^ta-qtc. 

When c = —a* and d = —b* (cf. (13.61) ). we set p = p* , q = q*, and 021 = <^Cii2 
so that Hn,m = cCnm- Moreover, if the constant coefficients of the three 
terms in Fn^rn are positive, then Fn^m > 0, so the complex conjugation re- 
duction is realized. Note that the above Fn,m can be written in a 2 x 2 
determinant form, 

{-(:;;::)( ( ^ ^ )( 
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Following the Gilson-Nimmo approach [42] in the continuous case, we 
construct the multidromion solution called the (M, A^)-dromion solution. 
The one-dromion solution corresponds to the simplest case of M = = 1. 
Hereinafter, we suppress the subscripts of the functions representing their 
dependence on the spatial variables n and m. When they are shifted, we 
express it using the shift operators 



■'n,m+l • 



We will consider various (M + A^) x (M + A^) matrices; they all have the 
same shape as 2 x 2 block matrices, so operations can be performed block- 
wise. For simplicity, off-diagonal zeros in the block diagonal matrices are 
omitted. We introduce two diagonal matrices as 



^2 



Si := diag(v3i, . . . , (pu), '^2 ■= diag(xi, ■■■,Xn) 
El := diag(^i, . . . , ^m), ^2 ■= diag(xi, ...,Xn) 



where 



We also introduce 

Q 

R :=\ ^ - 1, P := diag(pi, . . . ,pm), Q := diag(gi, . . . , giv) 



R :=[ ^ ^ ), P := diag(pi, . . . ,pAf), Q := diag(gi, . . . , gA^^ 



Q 
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Then, the following relations hold: 

»-'n" — "I Tit '-'m" — ^ 



I J ' "V g y ' 



^^^^-{^ g-i)^' ' O 

Note that the order of two diagonal matrices on the right-hand side can be 
changed because they commute. Moreover, we introduce M x M matrices 
K and Em as 

{K)ij :^ ^ ^ _ , {EM)ij:^l, l<i,j<M, 
1 - PiPj 

and N X N matrices L and Ej^ as 

{L)ki := ^ ^ _ , (£;iv)« 1, 1 < A), / < iV. 

They satisfy the relations, 

K - PKP ^Em, L- QLQ = En. 
We define (M + A'")-component column vectors as 

eM:= (1,..., 1,0,..., Of, e;v:=(0,...,0,l,...,lf, 

and 

I := Scm = ■ • • , (^M, 0, . . . , 0)^, m := Sejv = (0, . . . , 0, xi, • • • , Xjv)'^, 

7 := Scm = (^1, . . . , ^M, 0, . . . , 0)^, m := Sejv = (0, . . . , 0, xi, ■ ■ ■ , Xjv)^- 
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Then, we can easily show the following relations: 

Snl = Rl, Sml = I, S„m = m, Sm'n^ = R'n^, 
Snl = Rl, Sml = I, Snin = ifi, Sm'fri = R'fn- 
We set the "tau function" F^^^ as 

F = det J". 

Here, the (M + iV) x (M + A^) matrix J" is defined as 



where ^ is a constant (M + A^) x (M + A^) matrix. We also set the other 
"tau functions" Gn,m. and Hn,m as 

G = [ffi^J^-^Al) F, H = (T'^T-^Am}j F. 



Note that the one-dromion solution is reproduced by setting M = = 1, 
up to a minor redefinition of the parameters. Let us check that these "tau 
functions" indeed satisfy the only bilinear equation without time derivatives, 
fl3.17p . With the aid of the previous relations, we have 

= -AEeMeljE 
= -All^. 



Thus, we obtain 



Similarly, we obtain 



SnF = det (^J^-^ZF) 

= det (^J- J^-UzF) F 
= (^l-l^J^-'Al^ F. 

SmF = (1 - m^J^'^Am) F. 
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To compute SnSmF, we still need to know SnJ^ ^- This can be expressed 
as 

SnJ' = {SnJ') 



-1 



Here, we used the so-called Sherman-Morrison formula; recall that .4/ is 
a column vector and Z is a row vector. Combining the above results, we 
obtain 

{SnS^F) F={\-mF {SnJ'-') Am] {SnF) F 

= (SnF) (1 - fh^T-^Am) F - {fn^J^-^Al) F (^^J^-^Am}j F 
= {SnF) (SmF) - GH. 

This completes the proof of f l3.17p . 

It is a direct but lengthy calculation to check the remaining 12 equations 
fl3.18ap - fl3.21cp involving time derivatives. For example, in order to check 
fl3.18bp . we need the following intermediate formulas: 

dt^F = - {m^J^-^AR-'m) F, 

dt^H = (j^J^-^AR-'rnj {m^J^-^Am - l) F - (j'^^-^Arnj {rn^J^-^AR-'m) F, 
S-'H = (j^T-^AR-'m) F. 



To obtain the first formula, we first compute dt^J' as 



-AE 



O 



Q-'Ep 



-AR ^Scatc^S 
-AR-'mfnF . 



Then, we multiply it by J-" ^ and take the trace. 
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4 Concluding remarks 



In this paper, we have studied a suitable space discretization of the Davey- 
Stewartson system. The Davey-Stewartson system is an integrable NLS 
system in 2 + 1 dimensions, which involves two spatial variables on an equal 
footing and allows the complex conjugation reduction between the dependent 
variables. We started with a natural (2 + l)-dimensional generalization of the 
Ablowitz-Ladik lattice and then considered a nonlocal change of dependent 
variables to symmetrize the equations of motion. Consequently, we obtained 
the space-discrete Davey-Stewartson system inheriting most of the important 
properties of the continuous system; in particular, it is integrable and allows 
the complex conjugation reduction. The price to pay is the irrationality of 
the equations of motion and their high degree of nonlocality, which are not 
seen in the continuous case. Through a simple reduction, we reduced the 
degree of nonlocality and obtained a discrete modified KdV-type system in 
two spatial dimensions, namely, the (2 + l)-dimensional modified Volterra 
lattice fITTH]) . 

The (2 + l)-dimensional Ablowitz-Ladik lattice, as well as the space- 
discrete Davey-Stewartson system, is a superposition of four elementary 
flows that are mutually commutative. Naturally, the number of elementary 
flows is equal to the number of directions on the square lattice. Note also 
that both the (1 + l)-dimensional Ablowitz-Ladik lattice and the continuous 
Davey-Stewartson system can be written as a sum of two commuting flows. 
We conjecture that the (2 -\- l)-dimensional Ablowitz-Ladik lattice and the 
space-discrete Davey-Stewartson system possess four infinite sets of higher 
symmetries. As in the original Ablowitz-Ladik lattice [HI US], each set of 
symmetries could be generated from a single discrete-time system using the 
Maclaurin expansion in the step-size parameter. It would be interesting to 
provide a more precise description within the framework of the Sato theory, 
e.g., the discrete two-component KP hierarchy (cf. [19J). 

Using the Hirota bilinear method, we have constructed exact solutions 
such as the multidromion solutions of the (2 -|- l)-dimensional Ablowitz- 
Ladik lattice and the space-discrete Davey-Stewartson system concurrently. 
Their solutions can be obtained from the same set of bilinear equations, al- 
though their bilinearizing transformations are rather different (cf. fl3.7ap and 
fl3.13al) ). Note that fl3.13al) reflects the irrationality of the space-discrete 
Davey-Stewartson system that can, however, allow the complex conjugation 
reduction. The solutions were derived as the common solutions of the four 
elementary flows. On the level of the bilinear equations, the four flows look 
fully symmetric and stand on an equal footing. Thus, it is relatively easy to 
construct their common solutions despite the high number of bilinear equa- 
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tions. 
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